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The resistivity anisotropy unveiled in the study of detwinned single crystals of the undoped 122 
pnictides is here studied using the two-dimensional three-orbital Hubbard model in the mean-field 
approximation. Calculating the Drude weight in the x and y directions at zero temperature for a 
Q=(7r,0) magnetically ordered state, the conductance along the antiferromagnetic direction is shown 
to be larger than along the ferromagnetic direction. This effect is caused by the suppression of the 
d yz orbital at the Fermi surface, but additional insight based on the momentum dependence of the 
transitions induced by the current operator is provided, ft is shown that the effective suppression 
of the inter-orbital hopping d xy and d yz along the y direction is the main cause of the anisotropy. 
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Introduction. One of the most intriguing puzzles in the 
study of the Fe-based high temperature superconductors^ 
is the discovery of unexpected transport anisotropics in 
detwinned single crystals of doped and undoped AFe2Si~2 
(A = Ba, Sr, Ca)i2 Studies of the in-plane resistivity^ 
showed that the effect is the largest at low doping x~2- 
4% in Ba(Fei_ x Co a; )2As2, but it is present even in the 
undoped limit x=0 at low temperatures, i.e. in the mag- 
netically ordered state with wavevector (irfi). Recent 
studies^ for the undoped 122 materials have revealed a 
low-temperature anisotropy (defined as R = pb/p a — 1) 
R~0A, 0.35, and 0.09, for A = Ba, Sr, and Ca, respec- 
tively. This anisotropy is counter-intuitive because along 
the a-axis the spins order in an antiferromagnetic (AFM) 
arrangement, while along the &-axis they are ferromag- 
netic (FM). Intuition based on, e.g., double-exchange 
mechanisms for manganites would suggest that the FM 
direction should be less resistive than the AFM one. Op- 
tical conductivity measurements concluded that this un- 
expected anisotropy is caused by changes in the pop- 
ulations of the orbitals d xz and d yz at the Fermi sur- 
face (FS)^ in agreement with early mean- field studies 
where this unbalanced FS orbital population, without 
long-range orbital order, led to results compatible with 
photocmission techniques^ 

Several calculations have recently addressed the exper- 
imentally observed transport anisotropy. Using a five- 
orbital Hubbard model treated in a mean-field approxi- 
mation, and calculating the Drude weights via the Fermi 
velocities at the FS, results compatible with experiments 
were reported^ This agreement was observed in regimes 
where long-range orbital order is not present, and indeed 
the FS redistribution of spectral weight among the d xz 



ing on the three orbitals widely believed to be the most 
important in pnictides, namely d xz , d yz , and d xy , and 
also by identifying the electronic hopping amplitudes that 
cause the anisotropy. Our main results are that the ex- 
perimentally observed anisotropy clearly appears in the 
three-orbital model, in a state that is (7r,0) magnetically 
ordered, and it is mainly caused by the suppression of 
mier-hopping d xz -d yz processes along the FM direction. 

Models and methods. In this manuscript, the three- 
orbital Hubbard model for the pnictides at overall elec- 
tronic density n=4/3 (per site and per orbital) will be 
usediifii The hopping amplitudes that reproduce the FS in 
the paramagnetic state, with hole and electron pockets, 
were already provided and discussed in detail in Ref. Hoi . 
However, to help the readers in the understanding of our 
results, in Table U these intra- and inter-hopping ampli- 
tudes (in eV units) are provided again. From Table U 
note that the hoppings involving the d xy orbital, both 
intra-orbital and also inter-orbital with d xz and d yz , are 
the largest in value, inducing a large Fermi velocity in the 
regions of the FS where the d xy orbital dominates. This 
suggests that the d xy may play an important role in the 
anisotropy. The full three-orbital Hubbard model also 
contains a Coulombic on-site interaction (see Eq. (11) 
of Ref. [13), involving an intra-orbital Hubbard repulsion 
U, an inter-orbital repulsion U' , and a Hund coupling J, 
with the constraint U'=U-2J. The Hartree mean-field 
approximation used here has also been much discussed in 
previous literature and the reader is referred to Rcfs. 
H2I for details. The mean-field order parameters are the 
three electronic densities of each orbital, i.e. n XZl n 



and 



L xy j 



and the three magnetic moments m xz 
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and d yz orbitals caused by the (7r,0) magnetic order— is and m xy , and they are all determined via the minimiza- 



needed to understand the experimental results. Other 
calculations also for the five-orbital Hubbard model ar- 
rived to similar conclusions 

In this publication, the transport anisotropy found in 
experiments is revisited from the perspective of a sim- 
pler three-orbital Hubbard model.— Our goal is to refine 
the intuitive explanations given in Refs. 0@, by focus- 



tion of the Hartree mean-field energy. In the mean-field 
equations, the wavevector Q = (tt,0) is assumed. 

Let us focus now on the optical conductivity cr(w). Fol- 
lowing well-known computational studies of cr(w) in the 
context of the cuprates^ let us define first the param- 
agnetic current operators in the two directions as 
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TABLE I: Tight-binding (TB) hopping parameters of the 
three-orbital Hubbard model used in this manuscript. The 
energy unit is eV. The labeling convention is l=d mz , 2=d yz , 
3=d xy . The 13 and 23 hoppings are all affected by a 
factor (— 1)'*' = (-l) i *+ i » J with i=(i x ,i y ) being the la- 
bel of the Fe sites of a two-dimensional lattice. This 
modulation takes into account the two-Fe unit cell of the 
original FeAs layers^ The TB Hamiltonian is defined as 
HwT=^2 l i aPa t"f(4 iata ci+i,p,a + h.c.) ) where c\ a a creates an 
electron at orbital a of site i with spin projection a. i + 1 
denotes nearest and next-nearest neighbor sites to i. 



The optical conductivity and Drude weight in the y di- 
rection can be obtained and expressed similarly. N is the 
number of sites. In our calculation, the Dirac S functions 
are regularized as a Lorentzian i5(w)«(l/7r)e/(a; 2 + e 2 ) 
with a small but finite broadening parameter e. 

Results. One of the main results found in our study 
is shown in Fig. Q] where er(oj) in the two directions is 
shown for a state with magnetic order (tt,0). The val- 
ues of the couplings U and J are representative of the 
so-called "physical region" that was previously unveiled 
for the same three-orbital model^ In other words, by a 
comparison between neutron scattering and photoemis- 
sion experiments against mean- field results, in previous 
studies it was concluded that the three-orbital model has 
a "physical region" (where theory matches experiments) 
in the range i/~[0.7,1.3] and J/J7~[0.15,0.33]^i where 
the state is simultaneously magnetic and metallic as in 
pnictide parent compounds. Our a(tu) study is restricted 
to that "physical region" . 



3x= E E - U if( c t a ,a c i+i,^<r-h.c), 

(i,l=x,x+y,x-y) a,/3,a 

■iy= E E - U if ( C t a ,a C i+l ~ h - C -)> 

(i,l=y,x+y,-x+y) aJ3,a 

(1) 

while the kinetic energy operators are 

E J2 t n ( c L^ i+l ^ a + h.c), 

(i,l=x,x+y,x-y) a,/3,a 

f y= E E^( c U c m^ + ^-)- 

(i,l=y,x+y,-x+y) ujl.cr 

The total current, up to the first order term in the 
external field A.= (A x ,Ay), can be written as J x —{j x + 
f x A x )/N and J y =(j y + f y A y )/N. The real part of the 
optical conductivity in the x direction is given by>^ 

Rea xx (uj) = D x 6(uj) 

and from the <r(u>) sum-rule, it can be shown that the 
Drude weight in the x direction D x isH 

D* = (M-TM _ 1 V- l(0o|jx|^n)| 2 

2tt 2N E n -E ' (4) 

where <po is the many-body ground state (in this case the 
mean- field Q=(7r,0) state), and (f> n represents the many- 
body excited states, also produced in the mean-field cal- 
culation, with Eq and E n their corresponding energies. 








0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 
co(eV) 

FIG. 1: (color online) Example showing a(ui) in the "physical 
region"— of the three-orbital model (e=0.02). The unit of 
a(u) is e 2 /h. The couplings are U = 1.0 eV and J=U/4. The 
AFM direction (i.e. the x direction for magnetic wavevector 
(7r,0)) has a larger zero frequency conductivity than the FM 
direction, as in experiments. The FM direction also has a 
peak at a finite frequency ~J. 



Figure [T] shows that a(ui) for the three-orbital model 
is found to be in good qualitative agreement with exper- 
iments, namely at small frequency cj, where the Drude 
peak is located, the weight of this peak is larger along 
the AFM direction (the x direction) than along the FM 
direction. Similar results were obtained in the entire 
"physical region" , see Fig. [2] In addition, the finite fre- 
quency peak in the FM direction was found to scale with 
J. The ratio Dafm/-Dfm (i-e. D x /D y ) in the range of U 
shown in Fig.[2]varies approximately between 1.6 and 2.2, 
in qualitative agreement with results for the five-orbital 
modclii£ Thus, it is here concluded that the three-orbital 
models is sufficient to reproduce the d.c. conductivity 
anisotropy found in experiments^ 
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tion of an orbital weight redistribution at the FS is 
well established;^ with the d yz orbital suppressed for 
Q=(7r,0), it is desirable to develop a more intuitive un- 
derstanding of its influence on transport properties. For 
this purpose, the kinetic energy and current operators 
will be expressed in momentum space as: 



k a, {3, a k 

= E E 3?<V)<L, a , a <±,p,. = E^( k )' 

k (x,j3,<7 k 



(5) 



where a is the direction index (x, y). From Eq.(4), the k 
contribution (unfolded first Brillouin zone) to the Drude 
weight is defined as £> a (k)=Di ja (k)-Z)2,a(k), where 



DiM = (M - T a (k)\(bo) 
2tt 2N 

_D 2 ,a(k) _ f y^ fe (^o[ja(k)|0 n ) X (0 n |ja|0o) 
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FIG. 2: (color online) Drude weight/-7r vs. U in the "physical 
region" of the three-orbital model, at J=U /4. In this regime, 
the inequality Dafm>£ ) fm holds. As U increases toward the 
upper limit shown, the Drude weights in both directions are 
reduced due to increasing insulating tendencies.— 

For completeness, in Fig. |3][a) the population of the 
three orbitals is shown in the range of t/'s studied. From 
this figure, it is clear that there is no orbital order since 

the orbitals d xz and d yz are nearly identically populated. since summing over k leads to Eq. (4) . 
Further increasing U eventually leads to a regime of 
orbital-order^ but the opening of a gap renders the sys- 
tem insulating. It is important to note that Fig. [3ja) con- 
tains results obtained by integrating the orbital-selective 
density-of-states over all frequencies, while if the focus 
is only the vicinity of the FS, the orbital- weight redistri- 
bution phenomenon is observed.— As shown below, this 
redistribution is important to understand the anisotropy. 
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FIG. 3: (color online) (a) Charge density of each orbital vs. U 
in the "physical region" of the three-orbital model, at J=U/4, 
and for spin order (ir,0). (b) m a vs. U in the same U range, 
at J=U/A, and for spin order (n,0). 

In addition, Fig. GJb) shows that the magnetic moment 
in the range investigated is compatible with pnictides 
neutron experiments, ranging from ~0.25 Bohr magne- 
tons (/is) for the 1111 to ~1 hb for the 122 compounds^ 

Intuitive origin of the anisotropy. While the no- 



FIG. 4: (color online) D2,a(k)/-7r in the "physical region"— 
of the three-orbital model ([7=1.0 eV, J—U/4) and with spin 
order (n,0). For a discussion of the results see text. 

From Fig. |3|b) it is clear that the (ir, 0) pocket con- 
tributes significantly to D2. y since this wavevector region 
has sizable intensity. However, the contribution of the 
(0, 7r) pocket to D 2 , x is negligible, thus inducing the sig- 
nificant anisotropy observed in the overall Drude weight 
(note that the rest of the highly intense features in i?2,y 
arc 7r/2-rotated those of D 2 , x and thus do not contribute 
to the anisotropy). The reason is the orbital weight re- 
distribution at the FS: according to nesting scenarios, 
the (0,0) pocket (with mainly d xz and d yz character), 
interacts with the (jr, 0) pocket (mainly d xy and d yz ) 
when the magnetic wavevector is Q=(7r,0). This interac- 
tion needs to be intra-orbital^ hence the d yz states at 
the (it, 0) pocket are raised above the FS, while the d xy 
states are not. On the other hand, the (0,7r) states are 
not moved above the FS if Q=(7r,0), and hence this phe- 
nomenon does not occur for D 2 x . Note that the results 
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of Fig. [4]based on the Drude weights directly address the 
anisotropy in transport, and complements the analysis 
based on the orbital-weight redistribution.—^ Also note 
that a similar analysis of D\ y and D\ x (not shown) does 
not lead to the same clear anisotropy that i?2,a provides. 

To further simplify the understanding of the anisotropy 
evident in Fig. U let us now focus on the most relevant 
electronic hopping processes. Analyzing the values of 
the hopping amplitudes (Table 1), it is clear that those 
involving the d xy orbital (both inter- and intra-orbital) 
should be the most relevant since the other hoppings are 
much smaller in magnitude. The direct intra-orbital hop- 
ping d xy -d xy (tfi in Table 1) is not suppressed at the FS 
and should equally contribute to charge transport in both 
directions. Thus, the conductance should not drop to 
zero in any of the two directions due to this intra-orbital 
contribution. However, the inter-orbital hopping d xy -d yz 
is suppressed at the FS in a magnetic state (tv,0). This 
hopping occurs only along the y direction, as previously 
discussed^ On the other hand, the hopping d xy -d xz is 
not suppressed and can contribute to electronic hopping 
along the x direction. For these reasons, an asymmetry 
is expected between the x and y directions in transport, 
as found in Fig. Q] In addition, for k close to the (ir, 0) 
pocket the inter-orbital d xy -d yz hopping, which only ex- 
ists along the y direction, needs an excitation to con- 
tribute to <t(w) along the FM direction (y-axis) because 
d yz is suppressed at the FS. This observation justifies the 
presence of a peak scaling with J in ofm (Fig. 1). 

To transform the intuition developed above based on 
hopping amplitudes into actual transport properties, the 
Drude weight will also be decomposed according to the 
hoppings corresponding to the different orbitals, via the 
following definitions: 
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(Mjf + '£ a \M^n\ja\h) (7) 



Drude weight D^ a is obtained from Eq. (7) by replacing 
faf) + fp a by faa and j*P+jP<* by j% a . The several 
Drude components obtained by this procedure are in Ta- 
ble [III From this Table, it can be seen that the main 
anisotropy arises from the fact that D 43 ~0.125 is an or- 
der of magnitude larger than D^ 3 ~0.011, due to the FS 
suppression of the d yz orbital. The rest of the contribu- 
tions in Table [H] that are unrelated to the inter-orbital 
hopping involving d xy arc similar in both directions and 
are not relevant to understand the anisotropy. Actually, 
for the largest of those, the naive intuition suggesting a 
better conductance along the FM direction is satisfied 
since D y 3 >Dl 3 . 
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TABLE II: Drude weight/-^ decomposed into the different 
orbitals (l=d x ~, 2=d yz , 3=d xy ) of the three-orbital model 
working at [7=1.0 eV and J=0.25!7. Finding negative Drude 
weights in some cases is a well-known effect arising from 
differences of two large numbers in Eq. (4). 



Summary. A mean-field study of o~(oS) employing a 
three-orbital Hubbard model for the magnetically or- 
dered parent compounds of the pnictides has been here 
reported. In agreement with experiments, the conduc- 
tance along the AFM direction is shown to be larger than 
along the FM direction. The simplicity of this model al- 
lowed us to reduce this effect to an intuitive picture: (i) 
The AFM conductance behaves normally with a notori- 
ous suppression of its value as compared with the non- 
interacting limit due to spin scattering, in agreement with 
intuition, (ii) However, along the FM direction the dras- 
tic reduction in the weight of the d yz orbital at the FS 
leads to a large effective suppression of the d xy -d yz hop- 
ping and associated conductance along that FM direc- 
tion, causing the anisotropy found experimentally. 



for inter-orbital hopping (a^/3) (a=x,y). The operators 
in Eq. (7) arise from Eq. (2) via f a =J2 aiCe0 T^ and 
j a =^2 a a p 3a^- F° r the case of intra-orbital, the diagonal 
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